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ABSTRACT 

Radial velocities measured from near-infrared spectra are a potentially powerful tool to search for 
planets around cool stars and sub-stellar objects. However, no technique currently exists that yields 
near-infrared radial velocity precision comparable to that routinely obtained in the visible. We are 
carrying out a near-infrared radial velocity planet search program targeting a sample of the lowest- 
mass M dwarfs using the CRIRES instrument on the VLT. In this first paper in a planned series about 
the project, we describe a method for measuring high-precision relative radial velocities of these stars 
from iGband spectra. The method makes use of a glass cell filled with ammonia gas to calibrate the 
spectrograph response similar to the "iodine cell" technique that has been used very successfully in 
the visible. Stellar spectra are obtained through the ammonia cell and modeled as the product of a 
Doppler-shifted template spectrum of the object and a spectrum of the cell, convolved with a variable 
instrumental profile model. A complicating factor is that a significant number of telluric absorption 
lines are present in the spectral regions containing useful stellar and ammonia lines. The telluric 
lines are modeled simultaneously as well using spectrum synthesis with a time-resolved model of the 
atmosphere over the observatory. The free parameters in the complete model are the wavelength scale 
of the spectrum, the instrumental profile, adjustments to the water and methane abundances in the 
atmospheric model, telluric spectrum Doppler shift, and stellar Doppler shift. Tests of the method 
based on the analysis of hundreds of spectra obtained for late M dwarfs over six months demonstrate 
that precisions of ~5ms _1 are obtainable over long timescales, and precisions of better than 3ms -1 
can be obtained over timescales up to a week. The obtained precision is comparable to the predicted 
photon-limited errors, but primarily limited over long timescales by the imperfect modeling of the 
telluric lines. 

Subject headings: instrumentation: spectrographs — techniques: radial velocities — stars: low-mass 
- stars: individual: (GJ551, GJ699, GJ406, GJ442B) 



1. INTRODUCTION 

1.1. Motivation for NIR radial velocities 

The vast majority of known exoplanets have been de- 
tected with high-precision (i.e. a < 10 ms -1 ) time-series 
stellar radial velocity measurements. Invariably, these 
measurements have been made using spectra obtained in 
the visible wavelength region. This is natural for two 
reasons. First, arguably the most interesting targets for 
planet searches are solar-type stars (broadly defined here 
as dwarfs with F, G, and K spectral types), and these 
stars are brightest at wavelengths shorter than 1 fim. The 
spectra of solar-type stars at these wavelengths are also 
rich in deep and sharp spectral lines suitable for Doppler 
shift measurements. The second reason is that visible 
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wavelength spectrograph technology is much more ad- 
vanced relative to that required for spectrographs oper- 
ating in other wavelength regions. When fed by a mod- 
erately sized telescope, visible wavelength spectrographs 
using a cross-dispersed echelle design and CCD detectors 
can yield high- resolution and high signal-to-noise (S/N) 
spectra with a large wavelength coverage for a large num- 
ber of solar- type stars. The more advanced technology 
in the visible also includes well-established methods for 
high-precision calibration. The combination of currently 
obtainable data quality and calibration together make it 
possible to reach radial velocity precisions of ~lms _1 
for many solar- type stars, which is sufficient to detect 
planets down to a few Earth masses in short-period or- 
bits a round them (e.g. lHoward et aLll2009t iBouchv et al.l 
120091) . 

Despite the obvious utility of radial velocity measure- 
ments in the visible, measurements at other wavelengths 
are also desirable. The near-infrared (NIR, i.e. 1 - 
5 /im) in particular is an interesting wavelength region 
for radial velocity exoplanet studies. In contrast to 
solar-type stars, low-mass objects that are cooler than 
~ 4000 K are brightest at wavelengths of 1 /im and longer. 
Furthermore, stars cooler than ~ 3200 K are so faint 
at visible wavelengths that high-precision radial veloc- 
ity measurements are impossible for all but a few very 
nearby examples even with the best instruments on 8- 
10 m class telescopes. For example, the sample of 40 
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stars in a long-term radial velocity planet search pro- 
gram s pecifically targeting low-mass stars that utilized 
UVES (|Dekker et aLlf2000f) at the VLT only included two 
objects with masses be low 0.2 M^, and only th ree with 
masses below 0.35 M Q (jZechmeister et aT1l2009| ). 

Very low-mass stars are an interesting sample of poten- 
tial planet hosts notwithstanding their neglected status. 
They are useful for probing the correla tion between gas 
giant planet frequen cy and stellar mass (lEndl et al.l l2006: 
iJohnson et £012 007) , are the most numerous stars in the 
Gala xy (jReid et al.lbool iHenrv etalll2006t ICovev et al.l 
l2008f ). and exhibit a larger dynamical response to or- 
biting planets an d have closer-in hab itable zones than 
higher- mass stars (|Kasting et al.|[l993T) . The latter point 
means that much less precision is needed to detect po- 
tentially habitable planets around low-mass stars with 
the radial velocity method than is needed for solar-type 
stars. The ubiquity of mid to late M dwarfs and the 
closeness of their habitable zones suggests that the clos- 
est stars hosting tr ansiting habitable pl anets are likely 
of this stellar t ype |Deming et all [20 09). Furthermore, 
iDeming et al.l ( 20091 ) have estimated that basic atmo- 
spheric characterization of a few such planets could be 
obtained with the James Webb Space Telescope using 
the techniques of transit and occultation spectroscopy. 
Therefore, very low-mass stars are an important sam- 
ple for the possible future study of the atmosphere of a 
habitable planet. However, we first must be able to use 
the radial velocity method to identify or confirm (in the 
case the planet is first identified by a transit detection), 
and measure the masses of low-mass planets around very 
low-mass stars before these investigations can be carried 
out. 

As faint as the lowest-mass stars are at visible wave- 
lengths, their extreme redness means they are reason- 
ably bright at NIR wavelengths. For example, an M6 
star (M* ss 0.1 M Q ) at 10 pc will have a y-band magni- 
tude of ~ 15.5, which would make it a very challenging 
target for high-precision radial velocity measurements in 
this wavelength region with current telescopes and in- 
struments. However, the same star will have a -fT-band 
magnitude of ~8.5, which is bright enough that spectra 
with S /N suitable for high-precision radial velocity mea- 
surements could be easily obtained using existing high- 
resolution NIR spectrographs on 8-10 m class telescopes - 
especially those equipped with adaptive optics (AO) sys- 
tems. Therefore, if the issues and advantages of low- mass 
stars described above are to be studied and exploited us- 
ing the radial velocity method, then the NIR spectral 
region offers the only possible option without building 
larger telescopes. 

Radial velocities measured from NIR spectra also likely 
offer the advantage that activity induced "jitter" is re- 
duced relative to measurements in the visible. Magnetic 
activity plays a role in limiting the effectiveness of ra- 
dial velocity planet searches because it results in inho- 
mogeneities on the stellar surfaces (spots and plage). 
Visible spectra are thus modulated by stellar rotation 
and variations in the activity. These intrinsic variations 
mimic real radial veloci ty changes and are very confusing: 
for planet searches (e.g. iQueloz et al.1200 A iPaulson et al.l 
12004 iWrig ht 2005). It is thought that radial velocities 
measured from NIR spectra will exhibit reduced jitter 
relative to the visible because of the smaller contrast be- 



tween the cool inhomogeneities caused by activity and 
the average stellar surface at the longer wavelengths. Sig- 
nificant evidence exists to support the existence of this 
effect, and this evidence is discussed in the following sub- 
section. 

The issue of activity induced radial velocity jitter is 
particularly relevant for low-mass stars because it is well 
established that a much higher fraction of late-type stars 
are active than earlie r-type stars. In th e most compre- 
hensive survey so far, IWest et al.1 (|2004| ) found that the 
fraction of active stars increased from < 10% at spectral 
type M3, to 50% at M5, and ultimately 75% at the end 
of the main sequence. As high as these activity fractions 
are, they are likely only low er limits for nearby stars be- 
cause the IWest et al] (|200l ) survey included stars orbit- 
ing perpendicular to the Galactic plane, and such stars 
are potentially very old. Therefore, the fraction of active 
stars in a radial velocity pla net search of M d warfs will 
likely be even higher than the West et al.l (|2004l ) statistics 
suggest because they will be drawn from an intrinsically 
younger sample. This has serious implications for planet 
searches because it means that the majority of the lowest- 
mass stars will exhibit substantial activity induced jit- 
ter in radial velocities measured from visible wavelength 
spectra even if their faintness can be overcome. 

1.2. Previous work in the area 

The advantages offered by NIR radial velocities have 
been identified before a nd there has been so me previous 
work t o realize them. iMartfn et al.1 (|2006l ) used NIR- 
SPEC (|McLean et aIllr998F on Keck to make NIR ra- 
dial velocity measurements of the brown dwarf LP 944- 
20. Their measurements over six nights had an rms dis- 
persion of 360 ms -1 . This result was inconsistent with 
previously measured visible wavelength radial velocities, 
which exhibited a coherent periodic signal with an ampli- 
tude of 3.5kms -1 . This result is generally interpreted as 
confirmation that radial velocity jitter is indeed reduced 
at NIR wavelengths for active cool dwarfs. However, it 
should be noted that the NIR and visible data sets were 
not contemporaneous, and the possibility that the star 
was experiencing a period of reduced activity during the 
N IR measurements can not be ruled out. 

Zap atero Osorio et al.l (|2009f) have also reported ra- 
dial velocity measurements obtained from NIRSPEC 
data, in this case for the very low-mass star VB 10. 
Their data were obtained over seven years and have 
an estimated precision of ~ 300 ms -1 . These measure- 
ments seem to support the astrometr ic detection of a gi- 
ant pl anet around VB 10 claimed by iPravdo & z Shaklan 
(2009). However, our measurements of VB 10, which are 
presented in paper II and exhibi t a dispersion of only 
10ms " 1 , are inconsistent with the lZapatero Osorio et al] 
(2009) measurements and also appear to completely 
rule out the existe nce of the propo sed giant planet 
around this star dBean et al.l 120101) . In addition, 
lAnglada-Escude et al.l (2010) did not detect the expected 
reflex motion of VB 10 in their visible wavelength radial 
v elocities, which hav e a typical precision ~ 200 ms -1 . 

Bla ke et al] (|2007l ) report ed NIR radial velo city mea- 
surements with PHOENIX (jHinkle et al.ll2003D on Gem- 
ini South. They targeted L dwarfs and reached a pre - 
cision of 300 ms -1 over five nights. iPrato et al] ([2008) 
observed a sample of young stars and some late-type ra- 
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dial v elocity standard stars with CSHELL ([Greene et al.l 
[j~99"l on the IRTF for seven nights. They reached pre- 
cisions of ~ 100 ms -1 for the standard stars, and 100 - 
300 ms -1 for the young stars. The measurements for the 
young stars again indicated no coherence with significant 
periodicities seen in visible radial velocities, but al s o were 
not contemporaneous. Similarly, Huela mo et al.1 (j2008[ ) 
monitored th e young star TW Hya over six nights us- 
ing CRIRES ([Kaufi et al.H2004f ) on the VLT. Their mea- 
surements yielded an rms dispersion of 35ms -1 , which 
was inconsistent with a coherent signal having an ampli- 
tude of ~ 300 ms - 1 seen in visible rad ial velocities ob- 
tained pre viously dSetiawan et al.l I2008D and contempo- 
raneously. iSeifahrt fc Kaufll (|2008f) also used CRIRES to 
obtain continuous observations of an M giant over five 
hours and reached a precision of ~20ms . 

In addition to these stellar observations, there has 
been som e work on measuring NIR radial velocities for 
the Su n. iDeming et al.l (|1987t ) and iDeming fc Plvmatel 
( 1994) presented data from a long-term monitoring pro- 
gram using a Fourier Transform Spectrometer (FTS). 
Their i^-band measurements had a precisi o n of better 
than 5ms -1 . More recently. iRamsev et al.l ((2008.) used 
a laboratory grating spectrograph to measure the signa- 
ture of the Earth's rotational velocity in F-band spectra 
of integrated Sunlight. They maintained a precision of 
~ 10 ms -1 over a few hours in two separate experiments. 

1.3. Challenges of NIR radial velocities 

Despite the attention paid to obtaining NIR radial ve- 
locities of cool stars, no previous work has achieved a 
long-term precision on a star other than the Sun within 
an order of magnitude of the precision that is routinely 
obtained in the visible. One reason for this is simply that 
enough signal was not obtained to probe the true limit 
in some cases. However, bas e d on photo n statistics, the 
observ ations of iPrato et"aT1 ()2008l ) and iHuelamo et al.l 
(2008) should have at least achieved precisions around 
a factor of two better than they did. Another issue is 
that no studies with an external check on the obtained 
precision have been reported continuing for longer than 
seven days. Therefore, the obtainable long-term preci- 
sion in NIR radial velocities has not been systematically 
developed or evaluated. 

The main issue that has limited the pursuit of high- 
precision NIR radial velocities is the lack of a suitable 
calibration method. For example, u seful lines from avail- 
able emission lamps (e.g. ThAR, lLovis fc Pep"el l2007h 
are infrequent in the NIR relative to the visible. This 
is a problem for the existing high-resolution NIR spec- 
trographs because they generally have relatively short 
wavelength coverage for single exposures. Laser fre- 
quency combs have the potential to provide calibration 
lines in nearly any sp ectral region conceivable for stellar 
radial velocity work ([Murphy et al. 1l2007HLi et "alll2008t 
ISteinmetz et al.l l2008fT However, the technique is still 
not developed enough for long-term use at an astronom- 
ical observatory. Furthermore, none of the existing NIR 
spectrographs can be fed simultaneously with light from 
a calibration source like a lamp or laser comb system, and 
are not stabilized to the level necessary to do without si- 
multaneous calibration. In addition, there is also the 
issue that current NIR spectrographs experience signifi- 



cant variations in illumination on short timescales. This 
means that even if these instruments could be simultane- 
ously fed with light from a lamp, the calibration would 
not track all the necessary effects for high-precision ra- 
dial velocity measurements. Therefore, emission lamps 
and laser combs are not useful for high-precision calibra- 
tion in currently available NIR spectrographs, but could 
be for future instruments design ed for larger wavelen gth 
coverage and excellent stability (| Reiner s etaTll2010h . 

All of the previous studies discussed above that tar- 
geted stars other than the Sun utilized the telluric spec- 
trum imprinted on the data for in situ calibration, al- 
though the lines and methods util i zed va ried somewhat. 
In their study, Seifah rt fc Kaufi (2008) demonstrated 
that telluric N2O lines near 4.1 /jm were stable to a 
level of 10 ms -1 ov er five hours. On the other hand, 
IDeming et al.1 ()1987f ) showed that the telluric methane 
lines in the if-band varied with a semi-amplitude of 
20 ms -1 depending on the hour angle due to the pre- 
vailing winds over the observatory. In addition to the 
possible variability of the telluric lines, another issue 
with using telluric lines as a radial velocity fiducial is 
that measurements obtained over a full observing season 
could be affected by systematics because the stellar lines 
will move relative to the telluric lines due to the changing 
barycentric velocity of the observatory along the line of 
sight. This will lead to more or less blending of features, 
which could influence the measurements. 

Th e well known "iodine cell" method ([Butler et al.1 
11996( 1 has been used to measure radial velocities to preci- 
sions of a few ms -1 with numerous visible spectrographs 
that are not highly stabilized. The general gas cell tech- 
nique involves placing a cell in front of the spectrograph 
during observations so that a standard spectrum is im- 
printed on each obtained stellar spectrum. The cell's 
spectrum recorded during each exposure can then be 
used to precisely determine the spectrograph response at 
the exact moment of the observations and, thus, calibrate 
the simultaneously obtained stellar spectra for radial ve- 
locity measurements. 

The success of the iodine cell in the visible implies that 
a gas cell could be a useful way to calibrate currently 
available NIR spectrographs. The main constraint for 
the gas cell method is that the gas in question should 
exhibit numerous sharp and deep lines in a wavelength 
interval where the stars targ eted for radial velocity mea - 
surements have lines as well ([Campbell fc Wal ker 1979). 
Ideally, this region would also be free of atmospheric 
lines. However, this is a particularly challenging require- 
ment in the NIR due to the prevalence of strong atmo- 
spheric lines even in the traditional windows between wa- 
t er absorption b a nds. 

IDeming et ail (|1987ft and IDeming fc Plvmatel (|1994D 
used a cell filled with N2O to calibrate their FTS mea- 
surement s of integrated sunlight f rom 2.0 - 2.5 fim. In 
addition, Mah adevan fc Gel (j2009f ) considered a number 
of options for NIR gas cells and concluded that the gases 
H 13 C 14 N, 12 C 2 H 2 , 12 CO, and 13 CO together could pro- 
vide useful calibration in the 7i-band. However, there 
have been no reported uses of a gas cell for high-precision 
NIR radial velocity measurements for stars other than 
the Sun, and no such measurements obtained using a 
grating spectrograph. 
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We have recently developed a new gas cell for simulta- 
neous calibration of spectra obtained in the if-band, and 
are currently using it to conduct a radial velocity search 
for planets around very low-mass stars with the CRIRES 
instrument at the VLT. In this paper, which is the first 
in a planned series about our planet search, we describe 
a method for measuring high-precision relative radial ve- 
locities of cool stars using observations made with the 
cell, and we present extensive data sets demonstrating 
the performance of the method over long timescales. The 
paper is laid out as follows. In §2 we describe the new 
cell and its implementation in CRIRES. We describe the 
typical observing procedure and data reduction used in 
§3. We detail the radial velocity measurement algorithm 
in §4. In §5, we present tests of the obtained radial veloc- 
ity precision. We conclude in §6 with a discussion and 
outlook for future planet searches utilizing NIR radial 
velocities. 

2. AMMONIA CELL 

The cell that we have developed for measuring 
high-precision NIR radial velocities contains ammonia 
( NH3). Ammonia is a well established w avelength stan- 
dard for the NIR (e.g. lUrban et al.lll989l and references 
therein), but to our knowledge has not been used for 
calibration in astronomical observations. At room tem- 
perature ammonia is in its gaseous state and exhibits a 
rich molecular spectrum in the NIR even with the rela- 
tively low column densities possible in a cell to be used 
at an astronomical observatory. 

Although ammonia exhibits lines in different regions 
of the NIR, existing high-resolution NIR spectrographs 
yield relatively small wavelength coverage in a single ex- 
posure. Therefore, a careful consideration of the region 
for observations was necessary in parallel with the choice 
of a gas for the calibration cell. We chose to make obser- 
vations in a window in the A"-band because the very low- 
mass stars we are interested in exhibit numerous sharp 
and deep lines in this wavelength region suitable for ra- 
dial velocity measurements (mainly from the first over- 
tone transitions of CO), and ammonia exhibits a number 
of lines useful for calibration. The availability of a win- 
dow that can be spanned by existing NIR spectrographs 
and that contains a significant number of both calibra- 
tion lines and stellar lines is a particular advantage of 
ammonia over other possible gases, including those dis- 
cussed in §1.3. 

Unfortunately, the window we are observing in also 
contains a significant number of absorption lines arising 
from telluric methane and water. This was known be- 
fore we began our work. High-precision radial velocity 
measurements are usually made by avoiding regions con- 
taining telluric lines due to the expectation that these 
lines will exhibit variability on the order of a few to tens 
of ms ~ 1 . However, as discussed in §1, the lack of an ob- 
vious method for calibrating existing NIR spectrographs 
means that a more flexible approach is currently called 
for. We decided that using an ammonia cell in the K- 
band and accepting the contamination from the telluric 
lines was a good option considering all the competing 
issues. 

The main body of the ammonia cell we are using is 
a glass tube that has a length of 17 cm, and diameter 
of 5 cm. The windows were chemically fused on to the 



ends of the tube. They are made of quartz glass with 
a low OH content and have excellent visible and NIR 
transmission (the visible transmission is important for 
the use of an AO system) . Each window has a few degree 
wedge angle on one face to eliminate possible fringing 
when illuminated. The windows' wedges are oriented 
180° relative to each other to minimize the tilt induced 
on light passing through the cell. The cell was evacuated 
and filled with 50 mb ammonia at 15° C via a small tube 
in the side of the body. This tube was heated and pinched 
off while maintaining a constant pressure. In principle, 
the bonds holding the windows on to the body of the cell 
should remain sealed for 10+ years. 

We have implemented the ammonia cell in the 
high- resolution NIR spectrograph CRIRES (Kau fl et al.l 
120041 ). which is fed by the UT1 telescope of the Very 
Large Telescope facility. The cell is mounted inside the 
CRIRES "warm optics" box in an aluminum housing on 
a carriage that moves the cell in and out of the telescope 
beam. The cell sits just in front of the Nasmyth focus de- 
rotator in the converging f/15 beam from the telescope. 
At this location the cell is in front of all the spectrograph 
optics, as well as the instrument's integrated AO system. 
Observations of a star for radial velocity measurements 
are obtained with the cell in the beam, which causes the 
absorption lines of the cell to be imprinted on the stellar 
spectrum. The cell lines, whose position and shapes are 
well known, serve as a fiducial to precisely establish the 
wavelength scale and point spread function of the instru- 
ment at the time of the observation independently for 
each of the obtained spectra during analysis of the data 
(see §4). 

The ammonia cell is not temperature stabilized during 
use because this is not possible given the existing space 
limitations in CRIRES. Therefore, the cell experiences 
temperature variations driven by the ambient tempera- 
ture in the telescope dome during the observations. Its 
integration in the warm optics box provides a degree of 
thermal inertia, but does not prevent temperature vari- 
ations on timescales of hours or longer. A temperature 
sensor is used to record the temperature inside the box 
during every science exposure. The median temperature 
during our observations over six months has been 8.8° C, 
with a standard deviation of 1.7° C, minimum of 5.9° C 
and maximum of 11.9° C. Simulations using calculated 
synthetic spectra for ammonia suggest that a temper- 
ature difference of 10° C would result in a systematic 
radial velocity shift of ~ 1ms -1 . Our demonstration of 
5 ms -1 long-term precision with data taken using the 
cell (see §5) verifies that the temperature variations ex- 
perienced by the cell are at least smaller than this level. 
Therefore, the ammonia cell enables high-precision cali- 
bration without the need for additional heating or tem- 
perature stabilization equipment. Indeed, the fact that 
ammonia is a gas at the temperatures experienced at 
ground-based observatories, and that its spectrum is rela- 
tively insensitive to temperature variations seems to give 
it an advantage over other possible gasses for NIR cali- 
bration cells. 

We have made laboratory measurements of the cell us- 
ing the Lund Observatory Bruker IFS125 HR FTS to 
obtain the characteristic spectrum that is needed for the 
radial velocity measurements during analysis of the data. 
This FTS has a full optical path difference of 4.8 m, which 



NIR RADIAL VELOCITIES WITH AN AMMONIA CELL 



5 



TABLE 1 

Summary Of Observations 



Star exposure time (s) <S/N pixel 1 > exposures visit 1 # visits # epochs date range 

GJ551 20 160 12 27 4 2009.13 - 2009.61 

GJ699 20 150 12 12 4 2009.14 - 2009.61 

GJ406 60 140 8 11 3 2009.13 - 2009.46 

GJ442B 300 120 4 10 3 2009.13 - 2009.46 



gives a maximum resolving power of 10 6 at 5 /im. We 
used a resolution of Afc = 0.007 cm -1 (R = A/AA w 
620,000 at 2.3/xm), fully resolving the cell's spectral ab- 
sorption lines. The cell was illuminated with light from 
a 1,200.0 °C black body source. Broad band filters were 
used to select the wavelength region containing the am- 
monia lines. The wavelength scale was set by the HcNe 
laser used in the sampling of the interferogram, and 
the absolute accuracy for this region is on the order of 
0.001 cm -1 . The measurements were obtained with the 
cell cooled to 13° C. This temperature is close enough 
to the temperatures the cell experiences while in use 
in CRIRES given the stability of ammonia as described 
above (i.e. AT = 10° C corresponds to 1ms -1 ). The 
resulting spectrum has a measured S/N > 700 in the 
continuum near 2.3 /im. In addition to this spectrum, 
we also obtained measurements of the cell at 24° C on 
two separate occasions 13 months apart. Comparison of 
these data reveal no change in the cell's spectrum, which 
supports our assumption of the cell's stability. 

3. OBSERVATIONS 

We have been using CRIRES with the ammonia cell for 
an ongoing planet search targeting very low-mass stars 
since February 2009. All observations so far have been 
performed in four Visitor Mode runs with lengths be- 
tween four and six nights. We obtain observations of 
a portion of a single echelle order spanning the wave- 
length range 2292 - 2350 nm (a subset of the if-band), 
with three gaps of 3.0 - 3.5 nm that correspond to the 
spatial gaps between the four individual detectors in the 
mosaic used to record the data. The cell does not give 
enough suitable lines for calibration on the red-most de- 
tector. Therefore, only data from the three blue-most 
chips, which each contain enough ammonia lines, are 
utilized for the radial velocity measurements (total of 
36.4 nm spectral coverage). 

For most observations, including all those utilized for 
this paper, we use a spectrograph entrance slit with a 
projected width on the sky of 0'.'2, which gives a nomi- 
nal resolving power R rts 100,000. For the faintest stars 
we are observing (K > 9), a 0'.'4 slit is used to obtain 
more photon counts at the cost of reduced spectral res- 
olution. The CRIRES AO system is always utilized to 
improve the slit throughput. Multiple exposures (4 - 16) 
are taken over a time-span of up to 20 minutes during 
each of the observational visits to a star. The telescope 
is nodded along the slit during a visit so that half the ex- 
posures are obtained in one nod position, and the other 
half in another position. We limit the exposure times 
to <300s in order to minimize the possible variations in 
the sky background between data taken in different nod 
positions. The exposure times are set when possible so 
that the average S/N will be <~ ISOpixel" 1 . For the faint 



stars, the exposure time is simply set to its maximum 
value. An exposure time of 300 s typically yields S/N = 
lOOpixcl -1 for a late M dwarf having K=9.5. 

The obtained data are reduced and one dimensional 
spectra extracted using custom software that utilizes 
standard methods. Each chip is treated separately. First, 
the data are corrected for the signal arising from the bias 
voltage, dark current, and the amplifier glow by sub- 
tracting a dark frame taken with the same exposure time 
as the exposure in question. Multiple dark frames with 
the same exposure time are usually obtained during the 
course of an observing run. These are combined to create 
a master dark for the run using a sigma-clipped average, 
and this master dark is used for the correction. 

The next step in the reduction is the correction for the 
non-linear response of each individual pixel. The appro- 
priate corrections as a function of the measured counts 
were determined from analyses of sequences of flat field 
images taken with different exposure times. The correc- 
tions determined from data obtained at different epochs 
over the last year are consistent, which suggests that the 
detectors do not exhibit a significant time-varying re- 
sponse. 

After correcting for the non-linear response of the pix- 
els, the data are then corrected for the non-uniform re- 
sponse of the pixels. The obtained flat field images used 
for this are first corrected for the bias, dark, and ampli- 
fier glow signals, and non-linearity in the same manner 
as for the science exposures. The multiple flat field ex- 
posures taken during a run are then combined to give a 
master flat field frame using a sigma-clipped average. All 
the science data for the corresponding run are divided by 
this master flat. Bad pixels are flagged during the cre- 
ation of the master flat by identifying pixels whose sum 
over multiple exposures does not scale linearly with the 
exposure times. 

Once the data have been corrected for the detector ef- 
fects, the next steps are to subtract the background from 
the two dimensional images and extract one dimensional 
spectra. We do not utilize exposures obtained at different 
nod positions to estimate the background, although we 
have obtained data suitable for this procedure. Instead, 
the sky background for a given spectrum is estimated at 
each point in the dispersion direction using an average of 
the values of a number (20 to each side) of illuminated 
pixels that are more than l'/5 from the peak of the image 
of the star in the spatial direction (note that CRIRES is 
a long-slit spectrograph and is not cross-dispersed). At 
theses separations the count levels are always much less 
than 1% of the peak in the image of the star due to the 
use of the AO system to feed the spectrograph. We as- 
sume that the background level is the same for all pixels 
at the same point in the dispersion domain, which is a 
reasonable assumption despite the known curvature of 
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the slit because there are no sharp sky emission lines in 
the wavelength range we observe in. The background as 
a function of spectral position estimated from this pro- 
cess is then subtracted from the data to complete the 
reduction of the two dimensional images. 

We carried out tests of our method for estimating and 
subtracting the background compared to the method of 
simply subtracting an image taken in the opposing nod 
position. We found that a pixel-by-pixel subtraction of 
another science exposure leads to an increase in the noise 
consistent with that expected from the read noise of the 
detector. On the other hand, our method utilizes a num- 
ber of pixels (typically around 40) to assess the back- 
ground at each point in the dispersion domain, and thus 
introduces no significant noise while still delivering an 
accurate estimate. Therefore, this method yields slightly 
higher S/N in the end. 

One dimensional spectra and their corresponding un- 
certainties are extracted from the two dimensional re- 
duced data using an adaptation of the general o ptimal 
extraction algorithm described by iHornd ((198 6) . The 
background variance utilized for this algorithm is esti- 
mated as a function of spectral position during the cal- 
culation of the background level described above. This 
is usually slightly higher than expected from the known 
read noise of the detectors, but consistent with that esti- 
mated from propagating the uncertainties introduced at 
each step in the reduction. Pixels affected by cosmic ray 
strikes are identified and masked along with the previ- 
ously identified bad pixels during the spectral extraction. 
The final S/N for the spectra utilized in this paper is typ- 
ically 150 pixel -1 . A summary of the observations for the 
data presented in this paper is given in Table [TJ Included 
in the table are the nominal exposure times, typical S/N 
per pixel in the individual exposures, number of expo- 
sures obtained per visit, number of visits that have been 
performed, the number of epochs, and the time span of 
the observations. 

4. RADIAL VELOCITY MEASUREMENT METHOD 

4.1. Overview 

We measure the relative radial velocity of a star from 
a spectrum obtained with the ammonia cell using a 
me thod similar t o the "iodine cell" method pioneered 
by IButler et all (|1996f ). The basics of this method 
have been described extensively before (e . g . iValenti et al.l 
[19951: IButler "et~a^[l99fit lEndl et al.l[200l . Therefore, we 
only give a brief description of the overall method in this 
subsection, and the details of the unique aspects of our 
implementation in the following subsections. We refer 
the reader to these previous papers for more informa- 
tion on the theory and techniques of the general gas cell 
method. 

The core of the gas cell radial velocity measurement 
algorithm is the fitting of a model spectrum to an ob- 
served spectrum to simultaneously calibrate it and de- 
termine the Doppler shift of the stellar lines. The model 
is a composite of the different components that appear 
in the observed spectrum. In our case, spectra of the gas 
cell, star, and telluric absorption. 

The model for fitting the observed spectrum to deter- 
mine the stellar radial velocity is constructed by com- 
bining the different components in to a single spectrum, 
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Fig. 1. — Example model components and fit for the radial ve- 
locity measurements. The components of the model are given in 
the top three panels: the spectrum of the ammonia cell (top), syn- 
thetic telluric absorption spectrum (middle), and the template stel- 
lar spectrum (bottom). A comparison of a fit to the data is shown 
in the bottom two panels. (Top) The observed spectrum (points) 
and the best-fit model (line). (Bottom) The residuals from the fit 
(points). The short sections of spectra presented here represent 
~20% of the full wavelength coverage of one CRIRES detector. 

and then convolving this spectrum with an instrumen- 
tal profile (IP) and binning the result to the sampling 
of the observed spectrum. The free parameters in the 
model are adjusted to give the best match to the ob- 
served spectrum. The component spectra are combined 
in the model by sampling them on to the same wave- 
length scale and then multiplying them. Figure [T] shows 
an example of the model components and a fit to an ob- 
served spectrum. Our model calculations (multiplication 
and convolution) are done on the observed pixel scale 
oversampled by a factor of six. The critical aspect for 
the radial velocity measurement is that the stellar spec- 
trum is Doppler-shifted by an amount proportional to 
the radial velocity before the multiplication. This shift 
is one of the free parameters in the model, and its op- 
timization to give the best-fit to the observed spectrum 
constitutes the measurement. 

4.2. Template spectra 

Characteristic, or template, spectra of the three com- 
ponents that make up the model fitted to an observed 
spectrum during the radial velocity measurement pro- 
cess are required. The utilized spectrum of the gas cell 
is the FTS spectrum obtained at 13° C described in §2. 
The telluric spectrum used in the model is based on spec- 



NIR RADIAL VELOCITIES WITH AN AMMONIA CELL 



7 



tr um synthesis fo l lowin g exactly the methods described 
bv lSeifahrt etaLl <f2009h . In summary, the template tel- 
luric spectra are synthesized using a radiative transfer 
code with a time-resolved (3 hr resolution) model of the 
atmosphere over the observatory and an input line list. 
The radiative transfer code we use is LBLRTM_1 The in- 
put line list we use is based on the HITRAN database 
i|Rothman et al.ll2009fl . Initial tests using the HITRAN 
database as the input line list for the telluric spectrum 
synthesis resulted in poorer than expected fits to the ob- 
served spectra. This suggested that some of the HITRAN 
data could be modified so that the synthesized telluric 
spectra would more closely match the observed spectra, 
which would ultimately result in more precise radial ve- 
locity measurements. Therefore, we carried out work to 
determine improved data as described below. 

We determined improved line data for the telluric 
spectrum synthesis by fitting synthesized spectra to a 
CRIRES spectrum of a telluric standard star (a fast ro- 
tating early B-type star) that was obtained with the am- 
monia cell. The empirical spectrum exhibited no lines 
attributable to the observed star, and contained only tel- 
luric absorption lines and ammonia lines from the cell. 
The ammonia lines in this spectrum served as a fiducial 
for the determination of the new telluric line data in a 
similar fashion as they are used for the radial velocity 
measurements. We fitted the segments of the empirical 
spectrum from the three utilized CRIRES chips sepa- 
rately. That is, wc determined three unique line lists and 
imposed no constraints on the fitting linking the three 
lists to require, for example, the same methane or wa- 
ter abundance in the model atmosphere (see below), the 
same shift for lines that are part of the same band sys- 
tems, or etc. In the fitting for each chip, we determined 
an overall wavenumber shift for all the telluric lines in 
the spectrum. We also determined individual wavenum- 
bers and strengths for a total of 124 strong methane and 
water lines distributed fairly evenly over the three chips 
(i.e. the data for roughly 40 individual lines per chip 
were determined). 

The fitting process to determine the new telluric line 
data was similar to the process for the radial velocity 
measurement process, but without the stellar spectrum 
because the telluric standard star star exhibited no lines 
in the wavelength regions we observe in. For the fit, 
the abundances of telluric water and methane used in 
the model were allowed to vary to improve the match to 
the haze of weaker lines that could not be reliably an- 
alyzed individually. The wavelength scale and IP (see 
description in §4.3.1) were constrained by fitting the gas 
cell lines only. These lines were isolated by dividing the 
standard star spectrum taken with the gas cell by a spec- 
trum of the same star taken immediately before without 
the gas cell. The typical adjustment to the wavenum- 
bers of the telluric lines is on the order of 0.001 cm -1 , 
which is consistent with the estimated uncertainties in 
the HITRAN database and the uncertainty in the abso- 
lute scale of the gas cell FTS spectrum. The line data 
determination was done once, and the resulting line list 
is used for all the radial velocity measurements. 

The stellar template spectra we use for the radial ve- 

9 A description of the LBLRTM code and compilable source code 
are available at http://rtweb.aer.com/lblrtm_description.html 
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Fig. 2. — Example model IPs from each of the four observing 
runs when using a slit width of 0'.'2. 

locity measurements are based on CRIRES spectra of the 
same stars. We determine a template spectrum uniquely 
for each star. The template for a given star is created 
from an observed spectrum using spectral deconvolution 
in a manner similar to the standard practice in the iodine 
cell method. However, a key difference for our ammonia 
cell NIR work arises due the presence of the telluric lines. 

The process to create a stellar template begins with the 
obtaining of a spectrum of the star in question without 
the gas cell. Immediately before or after this observation 
we also obtain a spectrum of a telluric standard star at a 
similar airmass (typically within 0.02) and also without 
the cell. The standard star spectrum, which contains 
only telluric lines, is modeled to determine the IP and 
wavelength scale at the time of the observation. The stel- 
lar template spectrum is also divided by the the standard 
star spectrum to remove the telluric lines. The IP deter- 
mined from the analysis of the standard star spectrum 
is deconvolved from the stellar template spectrum after 
dividing out the telluric spectrum, and the determined 
wavelength scale from the analysis of the standard star 
is adopted. 

We do not utilize the gas cell in the process to create 
the stellar templates, which is in contrast to the typical 
technique for estimating the IP and wavelength scale of 
a template observation in the iodine cell method. The 
ammonia cell causes a different illumination of the spec- 
trograph than when it is not in use because of a small tilt 
in the beam induced by the wedged windows. Therefore, 
the spectrograph IP and wavelength scale is different for 
observations with and without the cell. The stellar tem- 
plate spectra are obtained without the cell because the 
goal of these observations is to get as pure of spectra of 
the stars as possible, and so in theory the best calibration 
exposures for these data would also be obtained without 
the cell. The spectra of the telluric standard stars taken 
without the cell contain only telluric absorption lines, 
and these lines offer the possibility of calibration because 
we have already developed an appropriate technique for 
modeling them. 

Tests indicate that using the telluric spectrum model- 
ing described above on its own yields reasonable IPs, and 
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a wavelength scale probably more closely matching that 
of the stellar template spectrum taken without the cell. 
Ultimately, the final radial velocities are insensitive to 
the choice of the exact method used for this step because 
the measured velocities are all relative and the absolute 
accuracy of the wavelength scale for the stellar template 
spectrum is not a limiting factor. 

4.3. Model parameters 

In addition to the stellar radial velocity, the other free 
parameters in the model are the wavelength scale and IP 
of the observed spectrum, continuum normalization, wa- 
ter and methane abundances in the model atmosphere 
used for the telluric spectrum synthesis, and Doppler 
shift of the telluric spectrum. 

4.3.1. Wavelength scale and instrumental profile 

The wavelength scale is parameterized by a second or- 
der polynomial. The IP model is constructed from the 
summation of 11 Gaussian pro files similar i n style to a 
model originally suggested by iValenti et al.l (119951) and 
described more recently by iBean et al.l (|2007T T In the 
IP model there is one central Gaussian, and five "satel- 
lite" Gaussians on both sides. The width of the central 
Gaussian and the heights of the satellites are the free 
parameters. The width and spacing of the satellites are 
fixed and set so that their half- widths just overlap. 

Some example IPs determined from the CRIRES data 
obtained with a 0'.'2 slit width are shown in FigureO The 
determined profiles are essentially Gaussians, but with 
some small deviations. The overall size and shape of the 
IPs also appears to vary with time. We believe both the 
deviations and time variability of the profiles are real. 
With regards to the first, the general morphology seen 
in the examples for different observing runs are consis- 
tent with the results from other spectra taken during the 
same run. The time variability is also observed when 
only a simple IP model consisting of a single Gaussian is 
used. Furthermore, the final measured radial velocities 
for known standard stars observed during our program 
(see §5) are significantly degraded when we use only the 
simple IP model for the analysis. Therefore, the devia- 
tions and time variability do not seem to be an artifact 
of using the complex, multi-parameter model. 

There are two reasonable explanations for the slightly 
non-Gaussian shape and time variability of the CRIRES 
IPs. First, using the AO system to boost the S/N often 
leads to under filling the entrance slit even when using 
a slit width of 0'.'2. This leads to asymmetries in the IP. 
The AO performance varies with the atmospheric condi- 
tions and, therefore, the IP varies with the conditions as 
well. A second issue is that there were significant varia- 
tions in CRIRES between the observing runs. For exam- 
ple, between the second and third run (May 2009) the 
instrument was warmed up and the detector mosaic was 
moved in the spatial direction so that it was centered on 
the slit. Also, a number of movable spectrograph com- 
ponents have poor position reproducibility (e.g. the en- 
trance and intermediate slits, grating, and prism). The 
set point for the temperature for various parts of the in- 
strument has been altered once as well. All these factors 
combine to suggest that indeed the CRIRES IP is slightly 
non-Gaussian and varying with time. The number of pa- 
rameters (i.e. the number of satellites used) in our IP 



model was set so that the model captures these effects 
without over-fitting the data. 

4.3.2. Telluric Doppler shift 

In addition to the standard model parameters in the 
gas cell method for determining high-precision radial ve- 
locities, we found that a Doppler shift must also be ap- 
plied to the telluric spectrum in our analyses. The pos- 
sible need for such a shift was anticipated because the 
stability of the telluric lines at the few ms _1 level was 
not known beforehand. However, the sizes of the shifts 
that came out of our analyses were unexpected. Over six 
months of observations we find that the analyses that 
yield constant radial velocities for known standard stars 
(see §5) also yield telluric spectrum velocities that vary 
by up to a few hundred ms _1 epoch to epoch (typical 
rms ~ 75ms -1 ). 

We carried out two tests to investigate the origin of the 
apparent telluric line velocity variations. The first test 
was an analysis of spectra for telluric standard stars that 
were obtained with the gas cell regularly during our pro- 
gram so far. These spectra do not have stellar lines and 
may be modeled using only the gas cell template spec- 
trum and a synthetic telluric spectrum. We analyzed the 
standard star spectra by carrying out the same kind of 
fitting as when determining the stellar radial velocity, but 
without that parameter. That is, we fitted the standard 
star spectra to solve for all the other model parameters 
described above except the stellar radial velocity. We 
analyzed 102 such spectra that were obtained over six 
months. We found that the obtained telluric velocity 
parameter was not varying at the same level as seen in 
the stellar radial velocity analysis, which is affected by 
the line blending. Instead, the determined telluric ve- 
locities have an rms of 15.8ms _1 . When the velocities 
were binned over each night, the rms is 7.3 ms -1 . These 
are interesting findings in their own right, and we plan to 
publish a more detailed study of the variability of the tel- 
luric lines in the future. For now it is sufficient to know 
that this preliminary analysis shows that for Paranal Ob- 
servatory, the telluric lines in our observing window are 
at least stable at a level of ~20ms _1 when they aren't 
confused with stellar lines. 

The second test we carried out to check whether the 
telluric line velocity variability was real or not was to 
determine the stellar radial velocities using a synthetic 
stellar spectrum (generated with the PHOENIX code) 
for the template instead of the one determined empir- 
ically as described in §4.2. The point of this test was 
that the synthetic stellar spectra do not have any possi- 
ble residuals from telluric line division or deconvolution 
as the empirical templates do. This test yielded similar 
results for the telluric velocity parameter as the first test. 

The results from the two tests described above sug- 
gest that the telluric lines are much more stable than 
indicated by our main analysis to determine the stellar 
radial velocities. It is clear that the large variability of 
the telluric lines is only seen when using the empirical 
stellar templates. Therefore, we believe that the large 
telluric velocity shifts seen epoch to epoch come from 
imperfect telluric line division when creating the empiri- 
cal templates. The imperfect telluric line division leaves 
residuals in the empirical stellar templates that are some- 
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Fig. 3. — Measured radial velocities for Proxima Centauri. The black points are the results for individual exposures, and the orange 
points are the data binned over a given visit (always 12 individual exposures except for the first visit of 16 exposures). The error bars 
for the individual exposures are not shown for clarity. Top left The results from a sequence of observations over two hours during one 
night. Another star was observed during the gap in the sequence. Top right The results from the first observing run. The blue dotted box 
indicates the data shown in the top left panel. Bottom The measured radial velocities obtained over six months (binned points only). The 
blue dotted box indicates the data shown in top right panel. 



what correlated with the telluric spectrum. These kinds 
of residuals are present despite our obtaining spectra of 
a standard star at a very similar airmass to the template 
observations for the telluric line removal. The residuals 
arise from the finite S/N of the observed spectra, and 
short timescale variations in the telluric absorption and 
the instrument. 

The telluric residuals in the stellar template produce 
an effect on the telluric velocity shift parameter via cross- 
talk with the real telluric lines in observed spectra that 
are modeled for the radial velocity measurements. Over 
an observing season, the stellar lines move relative to the 
telluric lines due to the changing barycentric velocity of 
the observatory along the line of sight to the star (tens of 
kms -1 peak-to-peak). Therefore, the stellar template in 
the model fit to the observed spectrum has to shift as well 
(this is what constitutes the measurement of the stellar 
radial velocity in the first pace). Therefore, the cross-talk 
between the template and observed spectrum requires 
that the telluric spectrum component in the model be 
shifted to absorb the correlated residuals during the fit- 
ting process. It is a subtle effect because the main drivers 
in the fitting are the match between the three component 
spectra in the model and their real corresponding lines 
in the observed spectrum. However, it is still necessary 
as either the telluric or stellar templates must be shifted 
to fit the data, and we have not been able to identify a 
way to filter out the effect in the fitting process. 

Tests of the obtained precision indicate that using the 
empirical stellar templates gives significantly better re- 
sults than synthetic templates despite the reduction in 
the degrees of freedom and the increased parameter un- 
certainties from correlations introduced by an additional 



free parameter (the telluric line shift). Despite being 
completely free of telluric residuals, the current state-of- 
the-art synthetic spectra for late M dwarfs do not suf- 
ficiently match the spectra of individual stars for high- 
precision radial velocity measurements. We are currently 
exploring the possibility of generating synthetic spectra 
specifically for this task, but it is uncertain whether such 
models will ever be more useful than the empirical tem- 
plates despite the latter's flaws. Currently, synthetic stel- 
lar spectra might be useful as a sanity check if we see 
radial velocity variations in the measurements made us- 
ing the empirical templates, but synthetic stellar spectra 
are not being used for the main velocity measurements. 
We have decided to use the empirical stellar templates 
and accept the telluric spectrum velocity parameter as 
a nuisance parameter for now. The issue of the telluric 
spectrum velocity parameter and its effect on the stellar 
radial velocities is taken up again in §5. 

4.4. Fitting the data 

The parameters that yield the best fit between the 
model and the observed spectrum being analyzed are 
estimated using a standard non-linear least squares al- 
gorithm. The adopted uncertainty in the determined ra- 
dial velocity is estimated from the covariance matrix and 
multiplied by the square root of the reduced chi-squared 
of the fit to account for the imperfect modeling. In addi- 
tion to this error estimate, we also calculate the intrinsic 
radial velocity content of the data by considering the 
characteristics of the stellar spectrum (i.e. how sharp, 
deep, and numerous are the spectral lines), the obtain- 
able level of calibration from the gas cell lines, and the 
S/N of the data. The intrinsic radial velocity content 
is the maximum level of precision that would be obtain- 
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able if the modeling process was perfect, and there was 
no covariance with other parameters - it is the photon- 
limite d precision. We e stimate this value using Eq. 6 
from iButler et al.1 ()1996l ) and calculate the obtainable 
precision for the stellar spectrum (i.e. the measurement 
precision) and the gas cell spectrum (i.e. the calibration 
precision). The results for the stellar spectrum and gas 
cell are added in quadrature to give an estimate of the in- 
trinsic radial velocity content of the data. Our estimated 
uncertainties in the radial velocities are typically a factor 
of two larger than the calculated intrinsic content of the 
data. 

The ammonia cell yields a spectrum similar in radial 
velocity quality to the spectrum of a slowly-rotating late 
M dwarf. Therefore, the calibration precision is on the 
same order as the Doppler shift measurement precision. 
Obviously this is not ideal. For example, iodine cells yield 
calibration precision a few orders of magnitude better 
than the intrinsic Doppler shift content of solar-type star 
spectra due to the very large numbers of lines molecular 
iodine exhibits in the utilized wavelength region. On the 
other hand, our situation is similar to that of using ThAr 
lamps for calibration in the visible where the calibrati on 
precision can be a limiting factor (jLovis fc Pepdl2007l ). 

The spectra from each detector during a single expo- 
sure are analyzed separately (they are not further subdi- 
vided) , and the final radial velocity for a given exposure is 
the weighted average of the results from the three detec- 
tors with usable data. The weightings are the estimated 
uncertainties from the fitting described above. The mea- 
sured velocities are corrected for the barycentric motion 
of the observatory along the line of sight at the mid- 
point of the exposures to give the final values that may 
be used to monitor the motion of the target. All of the 
stars we are monitoring are very close by, and many have 
a large proper motion. Therefore, some exhibit a signif- 
icant amount of secular acceleration. We calculate the 
ex pected trend for eac h star according the prescription 
of IKurster et al.l (|2003l) . and subtract it when necessary. 

5. VELOCITY PRECISION TESTS 

In addition to the very low-mass stars we are mon- 
itoring as part of our ongoing planet search, we have 
also frequently observed Proxima Centauri (GJ551) and 
Barnard's Star (GJ699) to test our observational and 
data analysis methodology. These stars are two of the 
few very low-mass stars for which it is possible to obtain 
high-precision radial velocity measurements in the visible 
wavelength range (due to their being very nearby, and 
thus bright even at visible wavelengths), and previous 
work has shown th em to be radial velocity constant at 
the level of 3 ms" 1 (IKurster et al. l l2003l lEndl fc Kiirsterl 
120081 : [Zcchmci ster et al.ll2009() . We also present here data 
from two of our planet search targets that so far appear 
constant. 

5.1. Results for Proxima Centauri 

Figure [3] shows the radial velocities measured for Prox- 
ima Centauri from spectra obtained during 17 nights over 
six months. These velocities are completely consistent 
with the previous findings that the star does not exhibit 
large radial velocity variations. The rms of the radial 
velocities measured for 100 individual spectra obtained 
over a six night observing campaign is 6.5ms -1 (see the 
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Fig. 4. — Residuals of the radial velocities for Proxima Centauri 
as a function of number of individual points binned. The circles are 
for the data taken during the first run (6 night span, 100 points 
total) and the solid line indicates the expected reduction in rms 
(N~ 1 ' 2 where N is the number of data points binned) normalized 
to the unbinned (N=l) rms. The squares are for all of the obtained 
data, and the dashed line corresponds to the expected reduction in 
rms for those data. 

top right panel in Figure 131 ■ The intrinsic radial veloc- 
ity content of the data (i.e. the photon-limited preci- 
sion) is typically 5.6ms" 1 , while our estimated errors 
are typically 9.1ms _1 . When the data are binned over 
complete visits (usually 12 individual exposures of 550 s 
total observing time) during this observing run, the ve- 
locities exhibit a dispersion of 2.6ms" 1 . The unbinned 
and binned velocities exhibit the same level of dispersion 
from timescales of two hours (see the top left panel in 
Figure [3]) up to six days. The expected dispersion in the 
binned data from the intrinsic radial velocity content is 
1.6ms -1 , while that expected from our estimated errors 
is 2.6ms" 1 . 

The rms of the individual data points over short 
timescales is 30% smaller than our estimated errors, but 
is similar to the intrinsic quality of the data (i.e. that 
expected from photon statistics). This suggests that our 
error estimation method might be slightly overestimat- 
ing the uncertainties on this timescale. The main differ- 
ence between our estimated errors and the photon limited 
precision is that we also take into account the imperfect 
modeling by adjusting the nominal errors according to 
the chi-squared of the fit to the observed spectrum. The 
main limiting factor in the modeling is the stellar tem- 
plate, which is not optimal because it is created through a 
deconvolution of a finite S /N spectrum. Spectral decon- 
volution in general amplifies noise and can not yield the 
exact underlying spectrum. Also, as discussed above, the 
empirical templates contain residuals from the telluric 
line division. Other secondary factors possibly leading 
to imperfect reproduction of the observed spectra are the 
limitations of the IP model and the synthesized telluric 
spectra. 

Over a short time-span the systematics arising from the 
imperfect modeling will be same because the IP, wave- 
length scale, telluric spectrum, and etc. will be similar 
for all the data. Therefore, all the measurements will be 
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Fig. 5. — Radial velocities for Proxima Centauri measured from the spectra on the blue-most CRIRES detector against some of the free 
parameters in the model used to fit the observed spectra. The colors of the points correspond to the different observing runs as indicated 
in the top left panel. 



similarly affected and the systematics will be reduced. 
However, the rms of the binned data over a few nights 
is slightly higher than expected from the photon limited 
precision, but consistent with our estimated errors. As 
the data are binned, and the expected precision drops to 
a few ms _1 , the imperfect modeling probably becomes 
an issue again. 

The dispersion in the radial velocities measured for 
Proxima Centauri over the full monitoring time-span (see 
the bottom panel in Figure[3]) is 11.7ms -1 for the indi- 
vidual spectra (not shown in Figure[3]) and 5.4ms" 1 for 
the binned data. The typical value for our estimated 
errors over the full time-span is 11.5 ms , and the ex- 
pected dispersion in the binned data from these errors is 
3.4ms" 1 . The typical photon-limited error is 7.0ms -1 , 
and the expected dispersion in the binned data from 
these errors is 2.0 ms -1 . 

The dispersion in the results for the individual spectra 
over long timescales is larger than that expected from the 
intrinsic radial velocity content of the data. However, the 
dispersion is consistent with our estimated errors and this 



suggests that the imperfect modeling is important over 
long timescales, which is consistent with the discussion 
above. However, the rms of the binned data over six 
months is 60% higher than expected from the estimated 
errors, which indicates that an additional source of noise 
is present on long-term timescales. We attribute this ad- 
ditional noise to the cross-talk between the telluric lines 
in the spectra and residuals from the telluric line divi- 
sion in the stellar template used to model the spectra in 
the determination of the radial velocities as described in 
§4.3.2. 

Figure [4] shows how the rms of the Proxima Centauri 
radial velocities depends on the number of points binned. 
Over short timescales (i.e. a single observing run, circles 
in the figure) the rms of the data nearly follows the ex- 
pected trend (i.e. the rms goes as N" 1 / 2 , where N is the 
number of data points binned - see the solid line in the 
figure). The noise floor over short timescales appears to 
be ~2ms -1 . It is worth noting that although Proxima 
Centauri is observable in the visible, the previously pub- 
lished data only indicate that it is at least stable at the 
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level of 3ms" 1 . It could be that the star is variable at a 
lower level, and 2 ms -1 is potentially only an upper limit 
to the noise floor of our method over short timescales. 

Over long timescales (i.e. an entire observing season), 
the rms of the Proxima Centauri radial velocities is a fac- 
tor of two larger than on short timescales, and the data 
show significant deviations from the expected reduction 
in rms with binning (squares and dashed line in Figure 

. We attribute this to the additional noise arising from 
the telluric cross-talk described above. The noise floor 
over long timescales appears to be ~4ms _1 from these 
data, which is consistent with the results seen for other 
stars (see below). 

In Figure [5] we plot the radial velocities for Proxima 
Centauri measured from the spectra on the blue-most 
CRIRES detector against some of the free parameters in 
the model used to fit the observed spectra. This figure 
demonstrates the typical range for the parameters, and 
allows a check for correlations. We find that the radial 
velocities show no correlation with the IP width, methane 
or water abundance in the model atmosphere used for the 
telluric spectrum synthesis (shown as the column density 
of these species), wavelength zero point, or temperature 
of the gas cell. Within a given run (epoch), the stellar 
radial velocities are correlated nearly one-to-one with the 
telluric spectrum velocities. This is expected because the 
variability of the telluric spectrum velocity parameter on 
the order of a few hundred ms _1 arises from changes in 
the stellar velocity shift on the order of a few kms -1 . 
Within a given run the barycentric velocity of the ob- 
servatory along the line of sight will vary by at most a 
few tens of ms' 1 during all the observations of a given 
star. The telluric residual cross-talk causes the telluric 
lines to be dragged with the stellar template spectrum 
when it is shifted in the model. Within a given run, the 
shift correlation is just essentially one-to-one because of 
the small scales involved. 

As can be seen in Figure [SJ the stellar radial velocities 
are not correlated with the telluric velocities over the four 
epochs that the data presented here were obtained. The 
key point is that freeing the telluric velocity parameter 
yields constant radial velocities for the known constant 
stars over long-term timescales. This validates our use 
of this parameter. However, the telluric lines in general 
limit the obtainable precision, and the need for this shift 
parameter is an aspect of that. 

Proxima Centauri (K=4.4) is very bright relative to 
the stars we are targeting for the planet search (K = 7 - 
10). Therefore, we were able to obtain many more spec- 
tra of it during each visit (typically 12 spectra in less 
than 10 minutes) than we are obtaining for most of the 
other stars. Also, we specifically obtained a large number 
of observations of Proxima Centauri to carry out the ve- 
locity precision tests presented here. To investigate what 
level of precision the Proxima Centauri data suggest for 
the stars in our program for which we are obtaining fewer 
spectra per visit, we used a simulation. 

We typically obtain four spectra per visit for the planet 
search targets with a similar S/N as the Proxima Cen- 
tauri spectra. For the simulation, we randomly drew and 
binned four radial velocity measurements from the Prox- 
ima Centauri data for each night of observation to create 
an example measurement sequence. The rms of each of 
the trial sequences was calculated to estimate the ob- 
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Fig. 6. — Histogram of the rms residuals for the simulation of 
precision obtainable with four radial velocity measurements per 
night for Proxima Centauri. The dashed line gives the median of 
the sample, which is 6.4ms" 1 . 
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Fig. 7. — Radial velocities measured for Barnard's Star (top), 
CN Leo (GJ406, middle), and GJ 442B (bottom). 



tained precision. We repeated the simulation 100 times 
and a histogram of the resulting rms values are plotted in 
Figure |5] The median rms of the 100 trials is 6.4ms -1 . 
This is likely a good estimate of the obtainable precision 
for a realistic set of observations on our planet search 
targets. 
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5.2. Results for Barnard's Star and others 

In Figure [7] we plot the radial velocities we have mea- 
sured for the late M dwarfs Barnard's Star, CN Leo 
(GJ406), and GJ442B over time spans of four to six 
months. As mentioned above, Barnard's Star is a known 
radial velocity constant star (minus its significant secular 
acceleration. | Kiirster et aflT 2003) at the level of 3ms" 1 . 
iReine rs (2009) presented radial velocities measured from 
visible spectra for CN Leo obtained over three consec- 
utive nights. These data were strongly affected by a 
large flare that occurred during the observations, but 
outside this event the radial velocities exhibit a disper- 
sion of < 10ms -1 . No high-precision radial velocities 
have been published for either CN Leo or G J 442B over 
long timescales. CN Leo and GJ 442B are typical targets 
in our planet search and we present data here only as 
an illustration of our obtainable precision. Future work 
will focus on detailed interpretation of the results in the 
context of the planet search. 

The radial velocities for all three stars shown in Figure 
[7] exhibit dispersions on the order of 5ms -1 . Eight to 
twelve measurements per visit were obtained and binned 
for Barnard's Star and CN Leo. We performed simu- 
lations using the data for these stars to investigate the 
precision obtainable with only four measurements per 
visit as we did for Proxima Centauri. In both cases we 
obtained typical dispersions of 4 - 6 ms -1 . 

The star GJ442B is very similar in characteristics to 
the main planet search targets we are following and the 
usual four measurements per visit were obtained and 
combined. This star has V=15.4 and, therefore, would 
be a very challenging target for obtaining high-precision 
radial velocities with existing visible wavelength spec- 
trographs. However, as an M5 dwarf it is very red and 
it has K=8.0. With these data we have demonstrated 
that by using CRIRES at the VLT with our new gas cell 
we can obtain measurement precisions of ~ 5 ms -1 in 20 
minutes total exposure time on such a star. 

6. DISCUSSION 

The general technique of measuring radial velocities 
from NIR spectra has recently been receiving a lot of 
attention as a possible important direction for the obser- 
vational s tudy of exoplanets for the reasons described in 
§1.1 fe.g. ILunine et afll2008f) . However, up to now the 
feasibility of obtaining NIR radial velocity precision sim- 
ilar to that routinely obtained in the visible was largely 
unknown due to uncertainty about a number of tech- 
nical issues. The main open question has been how to 
calibrate NIR spectra, with a second issue being whether 
the quality of NIR detectors is sufficient. 

We have described a new gas cell suitable for calibrat- 
ing spectra of cool dwarfs in the -ftT-band, and demon- 
strated long-term radial velocity precisions of ~5ms -1 
from data obtained with CRIRES at the VLT when using 
the cell. It is worth noting that the CRIRES detectors 
have severe cosmetic blemishes relative to CCDs. How- 
ever, the effects seem to be stable and we have been able 
to calibrate them out to an acceptable level. Further- 
more, the detectors were not even the highest-quality 
NIR detectors that were available at the time CRIRES 
was built (ca. 2002, |D orn et al1 l2004). and are certainly 
not competitive with the best detectors available as of 



this writing. Therefore, the issue of the quality of NIR 
detectors requires attention, but will not likely be a lim- 
iting factor for new instruments aiming at~lms -1 NIR 
radial velocity precision. 

Our ammonia gas cell has proven to be useful for NIR 
calibration down to a few ms" 1 precision despite not 
being temperature stabilized. Therefore, it could be use- 
ful for high-precision radial velocity measurements with 
other instruments not originally designed for excellent 
stability. The precision afforded by the cell and the mea- 
surement method we have described is already an order 
of magnitude better than the best previously obtained 
long-term NIR precision on a star other than the Sun. 
As a result, it now enables the possibility of searching 
for planets around a much larger number of very low- 
mass stars than was feasible before. For example, our 
planet search project includes 31 objects with estimated 
masses below 0.2 M Q , 22 of which have estimated masses 
below 0.15 Mq. A previous high-precision radial velocity 
planet search utilizing a visible wavelength spectrograph 
on an identical telescope (UVES on the UT2 telescope of 
the VLT) was only able to target two objects with esti- 
mated masses below 0.2 M , and onl y one of these two 
has a n estimated mass below 0.15 M Q (jZechmeister et akl 
[2001 . 

In addition to enabling the search for planets around 
more low-mass stars, the gas cell and radial velocity mea- 
surement algorithm we have developed also opens up a 
new frontier on the search for potentially habitable plan- 
ets. The orbital period range for a planet in the habitable 
zone around a star with a mass M = 0.10 M Q would be 
3-21 days (|Selsis et al.ll2007t) . The 5 ms" 1 precision 
obtainable with our method corresponds to the velocity 
semi-amplitudes induced by a 2.5 M ffi planet or a 4.6 M ffi 
planet on the inner and outer edges of the habitable zone 
respectively. Therefore, it should be possible to detect 
Super-Earth type planets in the habitable zones of very 
low-mass stars with a reasonable expenditure of observ- 
ing time on current facilities. 

Altogether our results have shown that obtaining high- 
precision NIR radial velocities is possible, and we see no 
technical reason why a level of precision of 1 ms -1 could 
not be achievable with a highly stabilized NIR grating 
spectrograph. However, the design of such a next genera- 
tion instrument will have to carefully weigh the availabil- 
ity of calibration against the availability of information 
(i.e. spectral lines) from the stars being targeted like we 
have done with our observing program. 

As discussed above, we are carrying out a planet search 
using CRIRES with the ammonia cell and this paper is 
the first in a planned series from the project. The next 
papers in the series will present high-precision NIR radial 
velocities of our M dwarf planet search targets. The data 
will allow us to detect planets around these stars with 
similar masses and orbital separations as those that can 
be detected around solar-type stars using state-of-the-art 
visible wavelength radial velocities. 
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